Multi-fluid potential in the loop cosmology 
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The scalar field can behave like a fiuid with equation of state — wp^, where w G [—1, 1]. In 
this Letter we derive a class of the scalar field potentials for which w = const. Scalar field with such 
a potential can mimic ordinary matter, radiation, cosmic strings, etc. We perform our calculations 
in the framework of the loop cosmology with holonomy corrections. We solve the model analytically 
for the whole parameter space. Subsequently, we perform similar consideration for the model with 
a phantom field {w < —1). We show that scalar field is monotonic function in both cases. This 
indicates that it can be treated as a well-defined internal time for these models. Moreover we 
perform preliminary studies of the scalar field perturbations with this potential. We indicate that 
non-Gaussian features are present admitting for the possible observational constraints of the model. 



I. INTRODUCTION 

In the framework of Loop Quantum Cosmology 
(LQC)[l|lone can introduce phenomenological Hamilto- 
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Here effects of the quantum holonomies have been intro- 
duced. The classical limit corresponds to the case p, ^ 0, 
what gives sin (/xc) //2 — > c. The parameter p, is in gen- 
eral function of the canonical variable p. In our consid- 
erations we choose p, — \J A/|p| where A — 2\f^-K^l\y It 
was shown that this leads to the proper classical limit and 
can be treated as a unique choice in the loop quantization 
class [1]. Based on the Hamilton equation / = {/, Tiphcn} 
we can derive equations of motion for the canonical vari- 
ables (c,p, TT^). The variables (c,p) can be related to 

the standard FRW variables (c, |p|) = {-ikvl'^ ,a^v'^l^). 
Subsequently, with use of the Hamiltonian constraint 
Wphon = we can derive modified Fricdmann equation 
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where k = SttG and 
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The equation of motion of the scalar field holds classical 
form 
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The special limit of the above theory is the case 
— 0, then Tiphcn ^ Wcff- It means that solution 
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of the equation traces the mean value (p) (for rigor- 
ous verification of the effective scenario in LQC see ReL 
0). It is however not necessarily the true for the general 
case with 1^(0) 7^ 0. However, even when we are not sure 
whether theory is truly effective one, it is interesting to 
investigate phenomenological consequences of the mod- 
els with potentials. Such studies have been performed for 
the broad class of potentials. In ReL [t'I model ([1]) with 
constant potential {y — K/8ttG) has been solved ana- 
lytically. Power law, cyclic and bicyclic potentials have 
been numerically studied in ReL Recently model 

with exponential potential has been studied, exhibiting 
sudden singularity behaviour . Other examples can be 
found in ReL [Mini- 

In this Letter we perform restriction for the scalar field 
energy density and pressure = wpcf, with w = const. 
In this situation scalar field can mimic any matter with 
equation of state p — wp where w = const and w £ 
[—1,1]. Based on this restriction we derive expression for 
the scalar field potential which fulfills this condition. We 
repeat these considerations for the case of phantom field 
{w < -1). 

Subsequently we highlight possible applications of the 
model. In particular wc consider scalar field perturba- 
tions and indicate its non-Gaussian features. This can be 
important from the point of observational constraints for 
the model. Namely since the obtained potential is non- 
quadratic, interaction between different modes is present, 
leading to the nonlinear effects in the CMB radiation. 
Therefore presented model can be potentially verified 
with astronomical observations. 

We show that scalar field with derived potential is 
monotonic function of coordinate time. This is impor- 
tant property since in quantum cosmological models we 
are looking for a well-defined intrinsic time. This is nec- 
essary since the Hamiltonian constraint H = leads to 
the time coordinate independent Wheeler-DeWitt equa- 
tion 7i|\E') — 0. Therefore wave function of the uni- 
verse 1 5') does not depend on time explicitly and some 
new intrinsic time must be introduced to trace quan- 
tum evolution. Monotonic scalar field is good candi- 
date. Namely performing canonical quantisation we re- 
place TTfj, TT^ = —idr),. Then matter part of the Hamil- 
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tonian is 
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Now (j) plays a role of time in analogy with the Klein- 
Gordon equation. 



II. MULTI-FLUID POTENTIAL 

Energy density and pressure of the homogeneous scalar 
field are expressed as follows 
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In this section we consider restriction — wp^ where 
w — const. This assumption leads to the expression 
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With use of equations (IS]),© and considered condition 
we obtain 



The above equation has solution in the form 



cosh 
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where — pc{l ~ w)/2. It is worth to mention that 
for (/) — ^ CX3 we recover V{(t)) oc cxp{— •\/247rG(l + w)(t)} 
which is well-known classical limit. In Fig. [l]we draw 
potential ITTI) for the different values of parameter w. 

Particular situation corresponds to the case = — 
{w = —1). Field can get such an equation of state only in 
the limit — > 0. Therefore this particular case must be 
excluded from considerations. Namely, we cannot mimic 
pure cosmological constant with time dependent scalar 
field. 

We must to stress that potential in the form 
1/ cosh^(a;) has already appeared in the different con- 
text in the loop cosmology. Namely, it was obtained 
by Singh [l^l for the scaling solutions in LQC dual to 
those in Randall-Sundrum cosmology. These considera- 
tions were performed for the model with a self-interacting 
scalar field and matter with a fixed equation of state w. 
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FIG. 1: Multi- fluid potential with the different values of the 
parameter w. In the figure the value 0o = was assumed. 



III. ANALYTICAL SOLUTIONS 

In this section we show analytical solution of the model 
with potential pdjl . For the considered case we can ad- 
ditionally derive equation 
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With use of continuity equation + 3H{p^ + 
(or equivalently equation (|4]) ) we directly obtain 



(12) 
= 
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Comparing this to the expression 



(13) 

with pT|) we obtain 



p(0) = Pc (cosh' [v/6^G(l + u;)(0 - 0o)] ) ■ (14) 

Similar solution, however not in such an explicit form, 
has already been derived in Ref. [12] ■ To obtain time 
t dependence we use expression = i±2il/((^) with 
potential (fTTj) to obtain 
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This leads to 

p{t) = Pc (1 + 6^Gpc(l + wf{t - to)^) ^ . (16) 

This solution can be also found in Ref. [ll|. In case of 
the w = — 1 solution has exponential form and can be 
found in the paper 

In Fig. [2] we show time dependence of the canonical 
variable p for different values of parameter w. Resulting 
dynamics is singularity free bounce. 

The time derivative of the field (p can be obtained from 
the equation ((32)l and has a form 
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FIG. 2; Evolution of the canonical variable p for the different 
values of parameter w. 



field energy density and pressure are expressed as follows 

P0 = + (19) 

P4> = -l<j>^-V {(!>). (20) 

From this we see that + < —(f)'^. The parameter 
w in the equation of state for the phantom field takes 
values w < —1. 

Phantom matter was introduced in loop quantum cos- 
mology in Ref. 16]. Subsequently phantom field in the 
loop quantum cosmology has been studied in Ref. [l^. 
Also the models with interacting phantom field have been 
investigated [11, [11, |23| . 

Based on the continuity equation we derive equation 
of motion 



From this we see that (f> is monotonic function of time 
(growing or decreasing). Choosing growing solution we 
have 



arcsinh [y/6nGpc{l + w){t — to)] 



We draw this function in Fig. [31 
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FIG. 3: Evolution of the scalar field (j> for the different values 
of parameter w. 

To complete above reasoning we verify obtained solu- 
tions in Appendix[Xl Namely we show there how to solve 
equations of motion without imposing the condition w — 
const in equation of state. Instead of this we solve equa- 
tions from scratch assuming the form of potential (fTT|) . 



IV. PHANTOM FIELD 
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Performing calculations in the same way like in the sec- 
tion |TT] we obtain 



cosh 
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where K = pc{l — w)/2. 

Solution for the scalar phantom field has a form 

m ^ 00 + "''^""^^ [^^6ggl + H(t - to)] ^ 
^-67rG(l + w) 

It is monotonically growing function too. Therefore 
phantom field with potential (I22p can be treated as inter- 
nal time. The canonical parameter p depends on value 
of the (h as follows 



p{4>) = Pc (^cosh^ ^-67rG(l + u;)(0 - (t)a) 
Rewriting equation (221) to the form 
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y/-67rG{l + w){(j) - 0o) = l+6T:Gpc{l+wf{t-tof 

(25) 

one can obtain explicit time dependence of the parameter 
p. We show this function in Fig. [H We see that, in 
contrast to the bouncing behaviour, here universe shrinks 
to zero for t — > ±oo. To investigate this state we consider 
time dependence of the energy density of the phantom 
field. We have 
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In this section we perform analo gou s consideration for 
the case of phantom field [H, [ij, Tisj . For this kind of 



what is shown in Fig. O We see that despite the variable 
p tends to zero energy density decreases too. 
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FIG. 4; Evolution of the canonical variable p for the different 
values of parameter w. 
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FIG. 5: Evolution of the energy density of the phantom field 
P0 for the different values of parameter w. 



V. PERTURBATIONS AND 
NON-GAUSSIANITY 

In the previous part of this Letter we have considered 
models with a homogeneous scalar field. Now we are 
going to take into account perturbations. We consider 
simplified model where scalar field is perturbed while 
background geometry holds its homogeneity. Moreover 
we aim rather to indicate further applications of the po- 
tential (fTTj) than to perform full analysis. Therefore con- 
siderations presented should be seen as an invitation to 
the more detailed studies. 

We split the scalar field 0(x, t) as 

0(x,t) =0(t) + <5</)(x,t). (27) 
Here 4>{t) is homogeneous mode defined as 
1 
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where Vq is some fixed fiducial volume. Then equation 
for the background part is 



dV{cj)) 
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Therefore solutions of equations of motion for the back- 
ground part are those derived in Section [IIII with = 0. 
First order perturbation of the full equation of motion 
leads to 



-m'^{t)S(t) = 



(30) 



which is valid for |(50/0| ^ 1. Here we have defined the 
mass term 
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(31) 



where 



B{t) = cosh^ ^&ttG{\ + w) 



l + QirGpcil + wft^. 

(32) 



Performing Fourier transform 



(27r)3/2 a{t) 



(33) 



and introducing conformal time drj = dt/a{t) we can 
write equation for the mode functions u(k, rj) as 



— u(k, rj) + + mlf,] u(k, 7?) = 0. 



Here we have defined effective mass 
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Based on the equation (|T4|) we can write solution for the 
scale factor as follows 

a(t) ^ac{l + 6TTGpc{l + w)h'^)^^ . (36) 
Applying this to the equation (|35p we obtain 
,1 



2 A-f^{5 + 2w -3w^) 



where 



3(l-hu;)2(8-h3u;-9w2) 
2 2 + 3w 
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The full analysis of the equation ([M)) with tti^j^ given 
above is behind the scope of this Letter. However we can 
immediately investigate its behaviour in the pre-bounce 
phase t ~* ~oo. In this limit 



lQ + 6w- 18u;2 1 



(39) 



where we have changed time for conformal. Advanced 
and normalised solution of the equation l|34p in consid- 
ered limit is 



(fe,,?) = y^^/^e'f (lH + ^)ijW(_^fc) 



(40) 
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where 



9 (9^2 - 2u; - 7) 
4 (1 + 3u;)2 ■ 



(41) 



This solution is valid for z^^ > or equivalently for 
w G [—1, —7/9] U {1}. Therefore only in the small range. 
Power spectrum of perturbations is then given as 
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where super-horizontal approximation 



ij(i)(^)^_lr(n)(|) " for x«l (43) 



has been used 
spectrum \ v 



It is worth to mention that scale invariant 

i_ I I is recovered only for w — —1. Therefore 

duality investigated in Ref. [2l| is broken here due to 
additional mass term. In the range w e (—7/9, 1) {v^ < 
0) we can obtain solutions by the replacement v = — ii>. 
Then for —kr]<^ 1 we have Vsc/) oc as discussed in Ref. 

Presented solutions are valid in the pre-bounce phase. 
Later when modes cross the bounce solutions of the equa- 
tion of motion have a different form. Despite the fact 
that size of horizon diverges at the turning point {H 
=> l/|ff| oo) the solutions evolve smoothly without 
any anomalies. In the context of loop cosmology such a 
behaviour has been recently approved for gravitational 
waves based on the analytical and numerical calculations 
[23i] . Interesting generic feature of the bouncing cosmolo- 
gies is suppression of power of the low energy models 
[2^ . [25} . Therefore we can expect damping of the low en- 
ergy modes also for the model studied here. However it 
can depend whether condition |(50/(/)| <^ 1 for validity of 
the linear approximation is conserved during the whole 
evolution. If this condition is not satisfied then nonlin- 
ear effects should be also taken into account. Otherwise 
suppression of power of the low energy modes should be 
easily obtained. Such an effect is typically presented as 
a possible explanation of the low CMB multipoles sup- 
pression [26j . However it seems that fine tuning must 
occur to match energy scales of the bounce with the cos- 
mological ones at which suppression of the CMB modes 
is observed. These issues have to be addressed in the 
further considerations. 

Scalar field with potential (|lip is non-Gaussian. 
Namely modes with different wave numbers interact one 
another Since observations suggest that primor- 

dial perturbations were nearly Gaussian [28] the CMB 
nonlinearities can in principle be used to restrict models 
with interacting potentials. In our case potential of inter- 
action has non-polynomial form therefore non-Gaussian 
features of the field should be significant. However since 
we do not know how many perturbations originate from 
the short quantum phase (where model studied here can 
be applied) we cannot infer about their impact on the 
CMB nonlinearities. In principle contribution from the 



phase of bounce can be affected by the succeed phase of 
inflation. Detailed studies of this issues are however be- 
yond the scope of this Letter. We devote for these issues 
another paper [2^ while now concentrate only on some 
preliminary estimations of nonlinear effects. 
The lowest order interacting term is given by 
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This term becomes important when dominates the mass 
one. Therefore ratio 



1 m 
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can be seen as a measure of non-Gaussianity. Here we 
have defined parameter 



f(t)=sgn(0\/-At(l 



^^^^J^. (47) 



2B(t)-3y B{t) 
We show this function in Fig. [51 At the points where 



-i.o -0.5 




FIG. 6: Parameter of non-Gaussianity ^ for w — 0. 

B{t) = I non-Gaussian part becomes dominant since the 
mass vanishes there. In the limits 



hm m 

r— >±oo 



w) 



(48) 



the parameter ^ stabilises. Above results indicate that 
non-Gaussian effects can be in fact produced. Especially 
interesting are points where B{t) — |, then nonlinear 
effects are dominant. These regions could be an effec- 
tive source of the cosmological non-Gaussianity. It is 
therefore in principle possible to constraint this model 
from CMB nonlinearities. This could bring observational 
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bounds on parameter pc or equivalently Barbero-Immirzi 
parameter 7 in this specific model. However to prove 
this supposition further detailed analysis is required. In 
principle bispectrum and resulting parameter of nonlin- 
earities /atl have to be determined. 



VI. SUMMARY 

In this Letter we have considered self-interacting scalar 
field theory in the universe with holonomy effects re- 
sulting from the Loop Quantum Cosmology. We have 
performed constraint for the scalar field in the form 
= wptfi where w =const. Based on this restriction 
we have derived expression for the scalar field potential 
which fulfills this condition. This condition can be ful- 
filled for w S (—1,1]. The case w = —1 is excluded. 
Such a state can be obtained only in the case and 
cannot be mimicked with (p ^ 0. We have solved con- 
sidered model analytically. Subsequently, we have per- 
formed analogous investigations for the phantom field 
w < —1. Potentials for scalar and phantom fields are 
related by the symmetry (1 -I- w) ^ — (1 + w). 

We have shown how barotropic matter can be intro- 
duced in LQC in the well-defined manner. In this ap- 
proach matter content is mimicked by the scalar field 
with derived potential (fTTj) . We have shown that field (f> 
is monotonic function of time in both considered cases. 
This result indicates that scalar field (as well as phantom 
field) with derived potentials can be moreover treated as 
intrinsic time. The models therefore seems to be good 
candidate to the further purely quantum considerations. 

We have performed preliminary studies of the pertur- 
bations of the scalar field with obtained multi-fluid po- 
tential. In the pre-bounce phase super-horizontal modes 
are typically damped. Namely the power spectrum is 
Vs^ oc k^. However to verify that this spectrum sur- 
vives during the bounce phase further numerical analysis 
should be performed. If super-horizontal spectrum does 
not change qualitatively this can lead to the suppression 
of the low CMB multipoles. We have also indicated that 
non-Gaussian perturbations are produced. Further anal- 
ysis of this issue may provide observational bounds on 
the model. Therefore beside the purely theoretical re- 
sults as the method of introducing barotropic mater in 
LQC also phenomenological consequences were found. 
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APPENDIX A: SOLUTION OF THE EQUATIONS 
OF MOTION - VERIFICATION 



In this appendix we are going to verify obtained so- 
lutions (|16ll8p solving from scratch equations of motion 
(|2l4p with potential (|lip . Without loose of generality we 
assume here (f>o — 0. Therefore we consider a theory with 
potential in the form 



cosh 



^67rG'(H-u;') 



(Al) 



where K = Pc(l — w')/2. Here we do not assume what 
is the interpretation of the parameter w' . It is just some 
fixed parameter and its meaning will be clarified later. 
Equation (|4|) with potential (jAip takes the form 



tanh(Q0) 
cosh^(a0) 



(A2) 



w here to sim plify notation we have defined a = 
\/ 67rG'(l -I- w'). To solve this equation we postulate so- 
lution in the form 



(\)(t) = — arcsinh (/3t) 



(A3) 



where /3 is some unknown constant. Applying this solu- 
tion, we rewrite equation (jA2[) to the form 



H = 



t 



1 + (3H-^ 



1 



(A4) 



Now let us calculate expression for the energy density 
©. With use of we obtain 



P - 



1 + I3H-^ 



0^ 
2a2 



(A5) 



Applying this to the modified Friedmann equation ^ we 
have 

^ . n [(3^ {2a^prt^ - l) + 2a\p, - V,)] {(3^ + 2V,a^) 
H = 



(A6) 



One can easily verify than equations (|A4[) and (jA6p are 
equivalent only if 



= 2a^{pc - K) = 6ttGpc{1 + w'f 
Therefore we obtain 

arcsinh [-/67rGp^(l + w')t\ 



v/67rG(l + w') 
Moreover equation (jA4[) can be written now as 

6TTGpc{l + w')^t 



1 



p 3{l + w')l + 6TrGpc{l + w')H^ 



(A7) 



(A8) 



(A9) 
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with solution Applying potential (|Aip together with solution (jASp we 

p(t) (l + GnGpcil + w' ft^) ^(^+^ . (AlO) 

Now one can interpret the parameter w' . Namely let us 
calculate expression for the coefficient 



w = w'. (A12) 

It is now transparent that solutions found in section IIIII 
are recovered. Thanks to relation (jA12p we have dw /dt — 



p i02 _ -j/^^) 

w = — = . (All) 0, therefore all found trajectories follow the w = const. 

P4> W + V{(t)) 
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